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Recent neutron scattering experiments have revealed that the generic form of the magnetic excitations in the 
high-Tc cuprates of wide range of doping has the so-called "hourglass" shape; it features both upward and 
downward excitations at the incommensurate (IC) momenta spanning from the resonance peak at the commen- 
surate momentum (ft, ft). We propose the two-component spin-fermion model as a minimal phenomenological 
model which has both local spins and itinerant fermions as independent degrees of freedom. Our calculations 
of the dynamic spin correlation function provide good agreement with experiments and show: (1) the upward 
dispersion branch of magnetic excitations is mostly due to the local spin excitations; (2) the downward disper- 
sion branch is from collective spin excitations of fermions; (3) the resonance mode is a mixture of both degrees 
of freedom. 

PACS numbers: 74.20,74.20-z,74.50 



The study of spin dynamics has been a key research interest 
since the discovery of the high-T c superconductor because it 
is expected that the spin correlation holds crucial information 
for the mechanism of the high-T c superconductivity (HTSC). 
For long the two main observations in the neutron scattering 
experiments of high-T c cuprates are (1) the incommensurate 
(IC) peaks at low energy or at quasielastic excitations^ and 
(2) the so-called resonance peak at commensurate wave vec- 
tor at relatively high energy (30 ~ 50 meV}4 In early ex- 
periments the IC peaks were observed only in the deeply un- 
derdoped lanthanum cuprate (LCO) and the resonance mode 
was the hallmark of the fully doped two layer yttrium cuprate 
(YBCO). However, later experiments reveal that both features 
appear in both cuprates, although sensitively depen ding on 
doping level. And with more recent experiemntP^^IlEl a 
unifying form of the magnetic excitations in the cuprate super- 
conductors has emerged as "hourglass" shape of excitations 
around the wave vector (1/2,1/2) (hereafter in units of 2%/a), 
in which the low energy IC excitations form the downward 
dispersion branch and the high energy IC excitations form the 
upward dispersion branch, and the two branches of excitations 
merge at the commensurate momentum and at the resonance 
frequency. 

It is a pressing question to understand the origin of this 
"hourglass" shape excitations. Theoretical proposals up to 
now can be classified into two groups: (1 ) theo ries based on 
the spin dynamics in the presence of stri pe a 10 ' n \ and (2) Fermi 
liquid type theories of itinerant fermions 12 13 . The key idea of 
the first group of theories is that the stripes formed by doping 
in the two dimensional Cu-O plane can split the commensu- 
rate spin wave excitations into two IC branches at the wave 
vectors (1/2 +£, 1/2) and (1/2 — e , 1/2 ) or at their symme- 
try rotated vector positions by x < — > y depending on the di- 
rections of the stripes. The dispersions from each branch of 
the two IC modulation cross at the commensurate wave vector 
(1/2,1/2) at a higher energy, which can be identified as the res- 
onance mode. This picture provides a qualitative explanation 
to the hourglass dispersion and the resonance mode. How- 
ever, this type of theories has difficulty to be extended to the 
higher doping regime where the presence and the nature of 



the stripes is questionable. The second group of proposals are 
itinerant fermion theories with interactiorPESl In this type 
of theories, the resonance mode and the downward dispersion 
can be obtained, but the upward dispersion branch is not yet 
satisfactorily reproduced. 

In this paper, we propose a two component spin-fermion 
model '"S as a minimal phenomenological model to provide a 
natural and unifying explanation for the above mentioned neu- 
tron experiments of high-T 6 . cuprates. In this phenomenolog- 
ical model, the minimal set of low energy degrees of free- 
dom are the spin wave excitations of local spins and the con- 
tinuum particle-hole excitations of fermions. A similar phe- 
nomenological theory is also known as one component spin- 
fermion model and has been intensively studied by Pines and 
coworkers.^! The key difference from the one component 
model, is the presence of the spin wave excitations directly 
from the local spins in addition to the usual collective spin 
density excitations from fermions. Despite a gross similarity 
to the one component model, we argue that the necessity of the 
local spins for the minimal model of HTSC is impelling from 
the neutron experiments. In a mixed momentum and real- 
space representation the corresponding Hamiltonian is written 
as 

H = £4(k)e(k)c«(k)+ £ gS(r).4(r)d apC p(r)+// s , 

k,a r,a,p 

(1) 

where the first term is the fermionic kinetic energy and the 
second term describes the coupling between local spins S(r) 
and the spin density of the conduction electrons. The last term 
represents an effective low-energy Hamiltonian for the local 
spins. When the local spins have a short range AFM correla- 
tion, the bare (before coupling to the fermions) spin correla- 
tion function has the general form as follows^. 

X i(q,£2)=x j(Q,0)-[l+^|q-Q| 2 -i2 2 /A| G ], (2) 

where Q the 2D AFM ordering vector, and the spin gap energy 
Asg and the magnetic correlation length cj combine to give the 
spin wave velocity v 4 = Asg ■ ^- Counting the coupling term 
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to one loop order, the dressed spin correlation functions of the 
model are written as follows. 

Xs\q,&) = 3Co^(q,a)-g 2 -Xo,/(q,a) < 3 > 
X7 1 (q»") = Xo}(q,^)-^ 2 -Xo,s(q,n) (4) 

where %o f is the noninteracting spin susceptibility of the con- 
duction band of the fermions and Xo.s is introduced in Eq.(2). 

Having two degrees of freedom in the model, two spin sus- 
ceptibilities Xs an d Xf should be calculated on equal footing. 
Previous studies of the local spin correlation embedded in the 
fermion bath 16 considered only the imaginary part of Xo.f ( 
Landau damping ) to damp the spin wave excitations in Eq.(3) 
and the real part of Xo.f is assumed either already included in 
the local spin dynamics described in Eq.(2)or its effect unim- 
portant. In fact, when the coupling g is weak, this approach is 
reasonable. But in the strong coupling limit (when the dimen- 
sionless coupling constant X = g 2 ■ Xo./(Q,0) • Xo,s(Q,0) ~ 
0(1)), including both the real and imaginary parts as in the 
above equations is crucial. In the strong coupling limit, both 
spin susceptibilities %s(q,i2) and %/-(q,il) become a mixture 
of both the local spins and the itinerant fermions and they as- 
similate to each other with increasing the coupling strength 
X. 

To make a contact with experiments, we use a tight binding 
model for the fermion dispersion 

e(k) = -2f(cos(£ x ) + cos(fc y )) -2f'cos(£ A ) -cos(^ v ) -/j. (5) 

For calculations in this paper, we chose t =-0.4 t, and /j = 
—0.8 If. The energy scale t and the choice of parameters t , fi 
will be discussed later with the numerical results. The fermion 
susceptibility Xo.f is calculated both in normal state (NS) and 
in superconducting state (SS) assuming a canonical d-wave 
pairing A(k) = Ao[cos(£ v ) — cos(fcy)]. 

Fig. 1(a) shows 7m%/(q,0) scanned along q = [h, 1/2) in 
the superconducting state. The superconducting gap Ao = 
0.2? and the bare spin gap A$g = 1-1/ is chosen; the phys- 
ical spin gap is strongly renormalized by coupling with 
fermions. The dimensionless coupling constant X = 0.80 is 
chosen. The main effect of the coupling is to pull down the 
bare spin gap A$g below the particle-hole excitation gap of 
Xo,/(q,£2) (~ 2Ao), which then forms a sharp resonance peak 
at Q = ( 1 /2, 1 /2). Centering from this resonance mode, both 
the downward dispersion branch and the upward dispersion 
branch span out. The origin of the upward dispersion is ap- 
parently from the local spin wave mode (see Eq.(2)) and the 
origin of the downward dispersion is the itinerant spin excita- 
tions of Xo.f- This fact is identified in Fig. 1(b) which shows 
7mXo,/(q,ii) scanned along q = (h, 1 /2) in the superconduct- 
ing state. The shape and strength of the downward whisker 
like excitations in ImXo.f is quite sensitive to the band struc- 
ture (controlled by t,t ', etc), Fermi surface curvature (con- 
trolled by fj, t ), and the size of the d-wave gap A(k). 

With the coupling strength X = 0.8, the dressed suscepti- 
bility X/(q,£2) obtains features of both the local spin suscep- 
tibility Xo.s an d the itinerant spin susceptibility Xo.f, an d the 
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FIG. 1: (a) The dressed spin susceptibility /m%y(q = (h, 1/2), £2) in 
superconducting state. Parameters are A^c = 1 . If , Ao = 0.2?, and X = 
0.8 (b) The bare spin susceptibility /m%o.y (q,f2) in superconducting 
state, (c) Im%f(q,Q,) in normal state (Ao = 0). (d) Im%Qj(q,£l) in 
normal state. 



behavior of the dressed susceptibility %s(q,i2) is qualitatively 
similar to the one of %/(q,£2) in this coupling strength. There- 
fore we show only the results of X/(q,ii) in Fig.l.; the exper- 
imental data of neutron scattering should be the contributions 
from both susceptibilities. With a smaller coupling strength 
(X < 0.5) the two spin susceptibilities retain more of their 
original characteristics of the local spin wave and the itiner- 
ant fermion susceptibility, respectively. 

Fig.l(c,d) are the same plots as in Fig.l(a,b) but in nor- 
mal state. First, the resonance peak becomes a completely 
overdamped mode having only a hump like structure in 
7mX/(q,£i). Second, the downward whisker like dispersion 
disappears because the free fermion susceptibility %o,/(q,^) 
in normal state has no such structure as seen in Fig. 1(d). 
Lastly, the upward dispersion remains almost similar to the 
case of the superconducting state. In experiments the up- 
ward dispersion in normal state appears more smeared^. This 
could be due to the effect of enhanced damping process at 
higher temperatures which is not included in our calculations. 
The results of Fig.l(a,c) successfully reproduce the main 
features of recen t neutron scattering experiments in high-T c 
cuprate J 4 ! 5 ' 6 ! 7 ' 8 ! 9 ! ie., the resonance mode in superconducting 
state, the hourglass shape of the upward and downward disper- 
sions, and their drastic change in normal and superconducting 
states. In particular, these results are in excellent agreement 
with the data of YBCO e #l. 

The mechanism of forming the resonance mode in our 
model is illustrated in Fig. 2. When the inverse of the dressed 
susceptibilities of Eq.(3) and Eq.(4) crosses zero (which oc- 
curs simultaneously in both susceptibilities), the dressed sus- 
ceptibilities develop a resonance mode, a bound state, or a 



3 



-Re t , „(0,£i) 

-g 2 Re x 0( (Q,n) 




(b) 


Re i' : ,..(Q,«) 

a 2 Re I or (Q.«) 

9 2 lm x„,(Q,") 









£2 (in unit of Tt 



FIG. 2: (a) Plots of bare susceptibilities Re% l s (Q,Q.), 

g 2 Re%oj(Q,Q.), and g 2 Im%Qj(Q,Q.), respectively in super- 
conducting state. Parameters are A$q = 1 . If , Ao = 0.2f , and X = 0.8. 
The vertical dashed line is a guide to the eyes indicating the position 
of pole at ft = 0.28f. (b) The same as (a) in normal state (ie. Aq = 0). 



completely overdamped mode depending on the presence and 
the strength of the imaginary part at the position of pole. 
In Fig.2 we plot separately Rex^^Q^Q.), Re Xoj(Q^)y and 
Im%o y(Q,£2) to make this point clear. Fig. 2(a) is the case of a 
superconducting state. By tuning the coupling constant g and 
the bare spin gap value A$g as chosen in Fig. 1(a), the pole 
°f X/.s(Q = Q:^) occurs at Q. res ~ 0.28f. At this frequency a 
small amount of Im%oj still exists just below the p-h excita- 
tion gap, making the pole a resonance mode. Fig. 2(b) shows 
the case of the normal state (Ao = 0) with the same parame- 
ters as in Fig. 1(c). The position of pole occurs at a little higher 
frequency (£2 ~ 0.5?) and there is substantial amount of damp- 
ing from Im%oj making an overdamped mode as shown in 
Fig. 1(c). With this mechanism of damping of the pole, even 
in the superconducting state, a small variation of the coupling 
or the superconducting gap, the intensity of the resonance 
mode can vastly change. Finally, we would like to make a 
remark about the unique character of the resonance mode in 
our model. The line of Re%Q l s (Q,Q.) in Fig.2 is not a simple 

inverse of a static potential (for example, jj in a RPA calcu- 
lation of Hubbard model) but it carries its own dynamics and 
spectral density. Therefore, the resonance mode formed by 
coupling of two dynamic susceptibilities %q t and Xo,s should 
carry the spectral density from both the local spin wave and 
the fermion particle-hole excitations. 

To make a comparison of our calculations with experi- 
ments, it is important to fix the energy scale of the model. The 
tight binding band of Eq.(5) is widely studied to fit the ARPES 
data and the estimate of t varies from 150 meV to 300 meV 
depending on the doping and different cuprate compounds^! 
Our calculation results are in good agreement with neutron 



experiments in terms of energy scale if we choose t ~ 150- 
180 meV. This value of t is somewhat low compared to the 
estimates from ARPES experiments. One possible reason of 
this discrepancy is that the extraction of t value from ARPES 
is carried by fitting the whole range of the quasiparticle dis- 
persion. As a result the high energy dispersion sets the overall 
energy scale and the value of t. However, for calculations 
of the spin susceptibility the high energy quasiparticle excita- 
tions are almost irrelevant and it is the low energy quasiparti- 
cle excitations near Fermi level which determine the structure 
of the low energy spin susceptibility. With this reasoning it 
is quite possible that the effective t value near FS is in fact 
much smaller due to a renormalization of a strong correlation 
effect. Also we use the chemical potential level /j = — 0.8f 
in our calculations. This value of /j corresponds to the band 
filling n = 0.503 in our model, which is even higher than the 
half filling (n = 0.5). In our phenomenological model the half 
filling has no special meaning because we do not include any 
strong coupling effect after the phenomenological Hamilto- 
nian Eq.(l). On the other hand the Fermi surface (FS) cur- 
vature is the single most important parameter to control the 
degree of incommensurability (IC) and the strength of the 
downward dispersion in the free fermion susceptibility %o,f 
as shown in Fig. 1(b). The curvature can be tuned with t and /j 
in our model. We decided to fix t = — 0.4f and tune the value 
of as a free parameter. Using a lower value of /j (smaller 
band filling), both the degree of the IC of the downward ex- 
citations and the p-h excitations gap edge have tendency to 
increase. Using a more realistic band with more tight binding 
parameters as t ,t , etc., it will be possible to choose more 
realistic band fillings. 

Fig.3.(a-d) show the constant energy scans of %f(q,Q.) in 
the superconducting state for Q. = 0.2f ,0.28f,0.6f, and 0.8f, 
respectively. Constant energy scans of neutron scattering data 
of YBCG 5 and LBCC 6 show peculiar patterns of IC peak po- 
sitions in (q x ,q y ) momentum space at different energy cuts. 
In particular, the 45 deg rotation of the patterns from a low 
energy scan (below the resonance energy) to a high energy 
scan draws special attention recently and several theoretical 
explanations have been proposed 11 ! i Results of Fig.3(a- 
d) demonstrate that the two component spin fermion model 
can consistently explain this phenomena, too. First, Fig. 3(c) 
is the scan of %/(q,i2) at the resonance energy, Q. = 0.28/ 
with the same parameters as in Fig. 1(a). It shows a very in- 
tense peak at (1/2,1/2) indicating a very sharp resonance not 
only in energy but also in momentum space. Fig3(a) and 
Fig. 3(b) are scans at higher energies than the resonance en- 
ergy and Fig. 3(d) is a scan of lower energy cut. We colored 
the highest intensity positions with black color to emphasize 
the clear patterns. The lower energy scan (Fig. 3(d)) shows the 
IC peaks at (1/2 ± 8, 1 /2) and (1 /2, 1/2 ± 8) forming a di- 
amond shape pattern. The higher energy scans (Fig3(a) and 
Fig3(b)) show that the IC peak positions at ( 1/2 ± 8, 1 /2 ± 8) 
and (1/2 ± 8, 1/2 =F 8) forming a square shape pattern which 
has the symmetry of the 45 deg rotated from the low energy 
pattern. The results of Fig. 3 excellently reproduce the ob- 
served patterns of the constant energy scan data of neutron 
experiments reported in YBCO^and LBCO^. 
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FIG. 3: Constant energy scans of Im%f(q,£l) at (a) Q, = O.St, (b) 
O. = 0.6f, (c) Q. = 0.28f, (d) Q, = 0.2f. Parameters are A,s G = hit, 
A = 0.2t, and X = 0.8. 



In our model we can trace the origins of the patterns. The 
low energy IC peaks and diamond shape pattern is basically 
a reflection of the band structure and d-wave superconducting 
gap. The high energy IC peaks and the square shape pattern 
has more complicated origin. At and above the resonance en- 
ergy the dressed spin correlation Xf is the result of a strong 



interplay between the local spin correlation and the itinerant 
spin correlation. Therefore the high energy scan pattern is de- 
termined by a subtle interplay/competition between Xo,s(q, ft) 
and %oj(q,Q.). The presence of IC peaks itself is the manifes- 
tation of the high energy spin wave dispersion spanning from 
the AFM wave vector Q; so the incommensurability increases 
with energy. However, whether the pattern becomes a square 
or a diamond shape has no universal mechanism. We tested 
various combinations of parameters t , pt, A^g, Ao and X. The 
low energy diamond shape pattern is quite robust within our 
model. As to the patterns of higher energy scans, although the 
square shape is the dominant one, it is not absolutely robust; 
with different parameters the diamond pattern often appears, 
too, at high energy scans. Therefore, we think that the 45 deg 
rotation of the IC peak patterns may not be a universal fea- 
ture of high-T c cuprates; it can change with doping and with 
different cuprate compounds. 

In conclusion, we proposed a two-component spin fermion 
model to explain the neutron scattering experiments in high- 
T e cuprates. With the two spin degrees of freedom of the local 
spins and the itinerant spins, our calculations of the spin cor- 
relations reproduced the essential features of the experiments: 
the hourglass dispersions, resonance mode, their changes in 
normal and superconducting states, and the IC peak patterns 
of the constant energy scans. With this success of the two- 
component spin fermion model to describe the neutron scat- 
tering experiments, the pressing question is now what the mi- 
croscopic theory is for the phenomenological two-component 
spin fermion model; specifically how the local spin degree 
of freedom survives after doping from the parent insulating 
cuprate compounds. 
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